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I summarize here the results of a global fit to the full data set corresponding to 535 days of data of the Super-Kamiokande experiment 
as well as to all other experiments in order to compare the most likely solutions to the atmospheric neutrino anomaly in terms of 
oscillations in the Ua — ► u e , v^—^Ut, and — > v s channels. 



oo 

On 

> 

o 

ON 



> 
o 



00 ■ 

a: 
i 

^ ■ 



X 



1 Introduction 



Atmospheric showers are initiated when primary cos- 
mic rays hit the Earth's atmosphere. Secondary mesons 
produced in this collision, mostly pions and kaons, de- 
cay and give rise tcuelectron and muon neutrino and 
anti-neutrinos fluxes El. There has been a long-standing 
anomaly between the predicted and,observed / v e ratio 
of the atmospheric neutrino fluxes 0. Although the abso- 
lute individual or v e fluxes are only known to within 
30% accuracy, different authors agree that the /v e 
ratio is accurate up to a 5% precision. In this resides 
our confidence on the atmospheric neutrino anomaly 
(ANA), now strengthened by the high statistics-sample 
collected at the Super-Kamiokande experiment □. The 
most likely solution of the ANA involves neutrino oscil- 
lations. In principle we can invoke various neutrino oscil- 
lation channels, involving the conversion of v^ into either 
v e or v T (active- active transitions p or the oscillatiomaf 
i/fj, into a sterile neutrino v s (active-sterile transitions )p'B. 
This last case is especially well-motivated theoretically, 
since it constitutes one of the simplest ways to reconcile 
the ANA with other puzzles in. the neutrino sector such 
as the solar neutrino problem □ as well as the LSND re- 
sult u and the possible need for a few eV mass neutrino 
as the hot dark matter in the Universe. 

The main aim of this talk^is to compare the — ► u e , 
v^i — > v T and the — ► v s transitions using the the 
new sample corresponding to 535 days of the Super- 
Kamiokande data. This analysis uses the latest improved 
calculations of the atmospheric neutrino fluxes as a func- 
tion of zenith angle, including the muon polarization ef- 
fect and taking into account a variable neutrino produc- 
tion point Ex 
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2 Atmospheric Neutrino Oscillation Probabili- 
ties 

The expected neutrino event number both in the absence 
and the presence of oscillations can be written as: 



N u = N uu + N„ 



where 



N, 



a/3 



da 



N e = N ee + N u 



(h,cos 9 V ,E V )P C 



(1) 



dE v d(cos Q v 
e(Ep)dE lJ dEpd(cos 8 v )dh 



a/3 



(2) 



for 



and Pap is the oscillation probability of vp — > v a 
given values of E v ,cos6 v and h, i.e., P a p = P(y a — > 
vp\E u , cos v ,h). In the case of no oscillations, the only 
non-zero elements are the diagonal ones, i.e. P a p = 5 a p. 

Here nt is the number of targets, T is the experi- 
ment's running time, E v is the neutrino energy and $ Q 
is the flux of atmospheric neutrinos of type a — fi,e; Ep 
is the final charged lepton energy and e(Ep) is the detec- 
tion efficiency for such charged lepton; a is the neutrino- 
nucleon interaction cross section, and 9 V is the angle be- 
tween the vertical direction and the incoming neutrinos 
(cos 6 U —1 corresponds to the down-coming neutrinos). In 
Eq. (^) , h is the slant distance from the production point 
to the sea level for a-type neutrinos with energy E v and 
zenith angle 9 U . Finally, the-slant distance distri- 

bution which is normalized to onetl. 

The neutrino fluxes, in particular in the sub-GeV 
range, depend on the solar activity. In order to take 
this fact into account in Eq. (|^), a linear combination 
of atmospheric neutrino fluxes $™ aa: and <i ) ™ m , which 
correspond to the most active Sun (solar maximum) and 
quiet Sun (solar minimum) respectively, is used. 

For definiteness we assume a two-flavor oscillation 
scenario, in which the v^, oscillates into another flavour 



either v L , 



The Schrodinger 



"fi ' "e j "/i ' "a " x 

evolution equation of the v^ — vx (where X = e, r or s 
sterile) system in the matter background for neutrinos is 
given by 



dt 



vx 



H. 



H^x 



A 1 
VX 
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H u =V a + 



Am 2 
4E V 



cos26» 



Am 2 

H x = Vx 777- cos . 



Hfj,x — —- 



4£„ 



sin 26* 



t ,x 



where 



v= V2G £ p 1_ Y 

Ve M (*e ^n) 

Here Gf is the Fermi constant, p is the matter density at 
the Earth, M is the nucleon mass, and Y e (Y n ) is the elec- 
tron (neutron) fraction. We define Am 2 = m 2 — m\ in 
such a way that if Am 2 > (Am 2 < 0) the neutrino with 
largest muon-like component is heavier (lighter) than the 
one with largest X-like component. For anti- neutrinos 
the signs of potentials Vx should be reversed. We have 
used the approximate analytic expression for the matter 
density profile in the Earth obtained in ref. E£l. In or- 
der to obtain the oscillation probabilities P a p we have 
made a numerical integration of the evolution equation. 
The probabilities for neutrinos and anti-neutrinos are dif- 
ferent because the reversal of sign of matter potential. 
Notice that for the — > v T case there is no matter 
effect while for the — > v s case we have two possi- 
bilities depending on the sign of Am 2 . For Am 2 > 
the matter efects enhance neutrino oscillations while de- 
press antineutrino oscillations, whereas for the other sign 
(Am 2 < 0) the opposite holds. The same occurs also for 
— > v e . Although in the latter case one can also have 
two possible signs, we have chosen the most usually as- 
sumed case where the muon neutrino is heavier than the 
electron neutrino, as it is theoretically more appealing. 
Notice also that, as seen later, the allowed region for 
this sign is larger than for the opposite, giving the most 
conservative scenario when comparing with the present 
limits from CHOOZ. 



3 Atmospheric Neutrino Data Fits 

Here I describe our fit method to determine the atmo- 
spheric oscillation parameters for the various possible 
oscillation channels, including matter effects for both 
— > v e and — > v s channels. The steps required 
in order to generate the allowed regions of oscillation pa- 
rameters were given in ref. □. I will comment only that 
when combining the results of the experiments we do not 
make use of the double ratio, R^/ e / 'Rffu, but instead we 
treat the e and /i-like data separately, taking into ac- 
count carefully the correlation of errors. It is well-known 



Table 1: Values of \ 2 an d confidence level for each experiment 
in the absence of oscillations. For unbinned data the number of 
degrees of freedom is 2, for combined binned data is 10 and for the 
full sample is 40. 10-CL represents the probability (in %) for no 
new physics. 



Experiment 


x 2 


CL (%) 


Frejus 


0.558 


24.3 


NUSEX 


0.39 


17.7 


1MB 


8.39 


98.5 


Soudan 2 


7.13 


97.2 


Kam sub-GeV 


12.7 


99.8 


Kam multi-GeV (unbinned) 


8.7 


98.7 


Kam multi-GeV (binned) 


18.7 


95.5 


SuperK sub-GeV (unbinned) 


23.2 


99.9991 


SupcrK sub-GeV (binned) 


32.2 


99.96 


SuperK multi-GeV (unbinned) 


9.88 


99.3 


SuperK multi-GeV (binned) 


22.9 


98.9 


All Combined 


102.2 


99.99997 



that the double ratio is not well suited from a statistical 
point of vieWpdu/3 to its non-Gaussian character. Thus, 
following ref. 



X 



i, J 



we define the \ as 

Nj h )-(a 2 da + al)j}-(Nj a -Nf), (4) 



where / and J stand for any combination of the ex- 
perimental data set and event-type considered, i.e, / = 
(A, a) and J = (B,j3) where, A,B stands for Frejus, 
Kamiokande sub-GeV, 1MB,... and a, (3 = e,fj,. In 
Eq. (||) Nj h is the predicted number of events calcu- 
lated from Eq. (0) whereas Nf a is the number of ob- 
served events. In Eq. (|j) <j 2 a and a 2 h are the error ma- 
trices containing the experimental and theoretical errors 
respectively. They can be written as 



o- a {A) p a/3 {A,B)ap{B) 1 



(5) 



where p a p(A : B) stands for the correlation between the 
a-like events in the A-type experiment and /3-like events 
in B-type experiment, whereas cr a (A) and ap(B) are the 
errors for the number of a and /3-like events in A and B 
experiments, respectively. |— . 

We compute p a p{A, B) as in ref. Ii3. A detailed dis- 
cussion of the errors^and correlations used in the analysis 
can be found in Reftl. We have conservatively ascribed a 
30% uncertainty to the absolute neutrino flux, in order 
to generously account for the spread of predictions in dif- 
ferent neutrino flux calculations. In Table |l| we show the 
values of x 2 an( i the confidence level in the absence of 
oscillation. Next we minimize the \ 2 function in Eq. ([|) 
and determine the allowed region in the sin 2 28 — Am 2 
plane, for a given confidence level, defined as, 



Xrain 



4.61(9.21) for 90 (99)% CL. (6) 
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Figure 1: y 2 . for fixed Am 2 versus Am 2 for each oscillation 
channel for Super-Kamiokande sub-GeV and multi-GeV data, and 
for the combined sample. Since the minimum is always obtained 
close to maximum mixing the curves for — * v s for both signs of 
Am 2 coincide. 



In Fig. [I] we plot the minimum x 2 (minimized with re- 
spect to sin 2 26) as a function of Am 2 . Notice that for 
large Am 2 £ 0.1 eV 2 , the x 2 is nearly constant. This 
happens because in this limit the contribution of the 
matter potential in Eq (||) can be neglected with respect 
to the Am 2 term, so that the matter effect disappears 
and moreover, the oscillation effect is averaged out. In 
fact one can see that in this range we obtain nearly the 
2 for the Vn — > v T and Vn — ► v s cases. For very 



same x 

small Am 2 ^ 10~ 4 eV 2 , the situation is opposite, namely 
the matter term dominates and we obtain a better fit 
v T channel, as can be seen by comparing 



for the — 

the — > v r curve of the Super-Kamiokande sub-GeV 
data (dotted curve in the left panel of Fig. |) with the 
— > v s and — > v e curves in the left panel of Fig. |l|). 
For extremely small Am 2 ?S 10~ 4 eV 2 , values x 2 IS quite 
large and approaches a constant, independent of oscil- 
lation channel, as in the no-oscillation case. Since the 
average energy of Super-Kamiokande multi-GeV data is 
higher than the sub-GeV one, we find that the limit- 
ing Am 2 value below which x 2 approaches a constant is 
higher, as seen in the middle panel. Finally, the right 
panel in Fig. [l] is obtained by combining sub and multi- 
GeV data. 

A last point worth commenting is that for the — > 
v T case in the sub-GeV sample there are two almost de- 
generate values of Am 2 for which x 2 attains a minimum 
while for the multi-GeV case there is just one minimum 
at 1.5 x 10~ 3 eV 2 . Finally in the third panel in Fig. [I] wc 
can see that by combining the Super-Kamiokande sub- 
GeV and multi-GeV data we have a unique minimum at 
1.6 x !CT 3 eV 2 . 
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Figure 2: Allowed regions of oscillation parameters for Super- 
Kamiokande for the different oscillation channels as labeled in the 
figure. In each panel, we show the allowed regions for the sub-GeV 
data at 90 (thick solid line) and 99 % CL (thin solid line) and the 
multi-GeV data at 90 (dashed line) and 99 % CL (dot-dashed line). 



4 Results for the Oscillation Parameters 

The results of our % 2 fit of the Super-Kamiokande sub- 
GeV and multi-GeV atmospheric neutrino data are given 
in Fig. ^. In this figure we give the allowed region of 
oscillation parameters at 90 and 99 % CL. One can notice 
that the matter effects are similar for the upper right 
and lower right panels because matter effects enhance 
the oscillations for neutrinos in both cases. In contrast, 
in the case of — > v s with Am 2 < the enhancement 
occurs only for anti-neutrinos while in this case the effect 
of matter suppresses the conversion in v^s. Since the 
yield of atmospheric neutrinos is bigger than that of anti- 
neutrinos, clearly the matter effect suppresses the overall 
conversion probability. Therefore we need in this case a 
larger value of the vacuum mixing angle, as can be seen 
by comparing the left and right lower panels in Fig. ||. 

Notice that in all channels where matter effects play 
a role the range of acceptable Am 2 is shifted towards 
larger values, when compared with the — * v r case. 
This follows from looking at the relation between mixing 
in vacuo and in matter. In fact, away from the resonance 
region, independently of the sign of the matter potential, 
there is a suppression of the mixing inside the Earth. As 
a result, there is a lower cut in the allowed Am 2 value, 
and it lies higher than what is obtained in the data fit 
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for the — > v T channel. 

It is also interesting to analyse the effect of combining 
the Super-Kamiokande sub-GeV and multi-GeV atmo- 
spheric neutrino data. Comparing the results obtained 
with 535 days given in the table above with those ob- 
tained with 325 days of Super-KamiokandcEl we see that 
the allowed region is relatively stable with respect to the 
increased statistics. However, in contrast to the case for 
325.8 days, now the — ► v T channel is as good as the 
— ► u ei when only the sub-GeV sample is included, with 
a clear Super-Kamiokande preference for the — » v T 
channel. As before, the combined sub-GeV and multi- 
GeV data prefers the — ► vx, where X = r or sterile, 



over the 



soluti' 



To conclude this section I now turn to the predicted 
zenith angle distributions for the various oscillation chan- 
nels. As an example we take the case of the Super- 
Kamiokande experiment and compare separately the sub- 
GeV and multi-GeV data with what is predicted in the 
case of no-oscillation (thick solid histogram) and in all os- 
cillation channels for the corresponding best fit points ob- 
tained for the combined sub and multi-GeV data analysis 
performed above (all other histograms). This is shown 
in Fig. | 



It is worthwhile to see why the za, 



channel is 



bad for the Super-Kamiokande multi-GeV data by look- 
ing at the upper right panel in Fig. |3[ Clearly the zenith 
distribution predicted in the no oscillation case is sym- 
metrical in the zenith angle very much in disagreement 
with the data. In the presence of — > v e oscillations 
the asymmetry in the distribution is much smaller than 
in the — ► v T or — > v s channels, as seen from the 
figure. Also since the best fit point for — > v s occurs 
at sin(20) = 1, the corresponding distributions are inde- 
pendent of the sign of Am 2 . 



5 Atmospheric versus Accelerator and Reactor 
Experiments 

I now turn to the comparison of the information obtained 
from the analysis of the atmospheric neutrino data pre- 
sented above with the results from reactor and accel- 
erator experiments as well as the sensitivities of future 
experiments. For this purpose I present the results ob- 
tained by combining all the experimental atmospheric 
neutrino data from various experimental. In Fig. ^ we 
show the combined information obtained from our anal- 
ysis of all atmospheric neutrino data involving vertex- 
contained events and compare it with the constraints 
from reactor experiments such as Krasnoyarsk, Bugey, 
and CH002lI, and the acceleratoij-experiments such as 
CDHSW, CHORUS, and NOMAdH We also include in 
the same figure the sensitivities that should be attained 
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Figure 3: Angular distribution for Super-Kamiokande electron-like 
and muon- like sub-GeV and multi-GeV events together with our 
prediction in the absence of oscillation (dot-dashed) as well as the 
prediction for the best fit point for Vp — > v s (solid line), — > 
Ve, (dashed line) and —> u T (dotted line) channels. The error 
displayed in the experimental points is only statistical. 



at the future long-baseline experiments now under dis- 
cussion. 

The first important point is that from the uppers 
right panel of Fig. [| one sees that the CHOOZ reactotl 
data already exclude completely the allowed region for 
the — > v e channel when all experiments are combined 
at 90% CL. The situation is different if only the com- 
bined sub-GeV and multi-GeV Super-Kamiokande are 
included. In such a case the region obtained is not com- 
pletely excluded by CHOOZ at 90% CL. Present accel- 
erator experiments are not very sensitive to low Am 2 
due to their short baseline. As a result, for all chan- 
nels other than v ^ — > v e the present limits on neutrino 
oscillation parameters from CDHSW, CHORUS and NO- 
MAD Ej are fully consistent with the region indicated by 
the atmospheric neutrino analysis. Future long baseline 
(LBL) experiments have been advocated as a way to in- 
dependently check the ANA. Using different tests such 
Long-baseline experiments now planned at KEK_(K2K) 
Fermilab (MINOS)@ and CERN ( ICARUSB, NOE 
and OPERAE3) would test the pattern of neutrino os- 
cillations well beyond the reach of present experiments. 
These tests are the following: r appearance searches, 

{NC/CC) near and the 



NC/CC ratio which measures (j\f(jj(j(j^ 
muon disappearance or CC near /CCf a r test. 
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Figure 4: Allowed oscillation parameters for all experiments combined at 90 (thick solid line) and 99 % CL (thin solid line) for each oscillation 
channel as labeled in the figure. We also display the expected sensitivity of the present accelerator and reactor experiments as well as to 
future long-baseline experiments in each channel. The best fit point is marked with a star. 
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The second test can potentially discriminate between 
the active and sterile channels, i.e. — > v r and — > v s . 
However it cannot discriminate between — > v s and 
the no-oscillation hypothesis. In contrast, the last test 
can probe the oscillation hypothesis itself. Notice that 
the sensitivity curves corresponding to the disappearance 
test labelled as KEK-SK Disappearance at the lower 
panels of Fig. || are the same for the — > v T and the 
sterile channel since the average energy of KEK-SK is 
too low to produce a tau-lepton in the far detector. 

In contrast the MINOS experiment has a higher av- 
erage initial neutrino energy and it can see the tau's. 
Although in this case the exclusion curves corresponding 
to the disappearance test are in principle different for 
the different oscillation channels, in practice, however, 
the sensitivity plot is dominated by the systematic error. 
As a result discriminating between — > v r and — ► v a 
would be unlikely with the Disappearance test. 

In summary we find that the regions of oscilla- 
tion parameters obtained from the analysis of the atmo- 
spheric neutrino data on vertex-contained events cannot 
be fully tested by the LBL experiments, when the Super- 
Kamiokande data are included in the fit for the — > v T 
channel as can be seen clearly from the upper-left panel 
of Fig. [|. One might expect that, due to the upward 
shift of the Am 2 indicated by the fit for the sterile case, 
it would be possible to completely cover the correspond- 
ing region of oscillation parameters. This is the case for 
the MINOS disappearance test. But in general since only 
the disappearance test can discriminate against the no- 
oscillation hypothesis, and this test is intrinsically weaker 
due to systematics, we find that also for the sterile case 
most of the LBL experiments can not completely probe 
the region of oscillation parameters indicated by the at- 
mospheric neutrino analysis. This is so irrespective of 
the sign of Am 2 : the lower-left panel in Fig. [| shows the 
— > v s channel with Am 2 < while the — > v s case 
with Am 2 > is shown in the lower-right panel. 
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